Evolution speed in some coupled-spin models 
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We investigate the time evolution of some models with N spins and pairwise couplings, for the case 
of large N, in order to compare evolution times with "speed limit" minima derived in the literature. 
Both in a (symmetric) case with couplings of the same strength between each pair and in a case of 
broken symmetry, the times necessary for evolution to a state in which the simplest initial state has 
evolved into a nearly orthogonal state are proportional to 1/N, as is the speed limit time. However 
the coefficient in the broken symmetry case comes much closer to the speed limit value. Introducing 
a different criterion for evolution speed, based on macroscopic changes in occupation, we find a 
corresponding enhancement in rates in the asymmetric case as compared to the symmetric case. 

PACS numbers: 03.65.-w, 03.65.Ud, 03.67.-a 



INTRODUCTION 

In a recent series of papers, Giovannetti, Lloyd, and 
Maccone Q-Q have elucidated the role of entanglement 
in approaching theoretical minima of the times required 
for quantum systems to progress from a prescribed initial 
state to a state that is orthogonal or nearly orthogonal 
to the initial state. The criterion for evaluating the out- 
come is the reduction of the square of the overlap of the 
initial and final wave-functions to some value less than e, 
where e is fairly small, but not necessarily infinitesimal, 
P(t) = (*(t)|*(0)) < e. The "theoretical minimum" of 
evolution times that we refer to depends on the lesser of 
the two quantities, 



AE = V(# 2 )-<#) 2 



(1) 



and E — Eq =< H > —Eq, where Eq is the ground state 
energy of the system and the brackets <> in all cases 



stand for the expectation value in the initial state. In 
terms of these quantities the fastest evolution time for 
meeting the criterion P(t) = e, as derived in refs. H-0 
and sharpened in refs. Q- Q is given by, 1 



t > max 



7ra(e) Tr/3(e) 



2(E-E ) 2AE 



(2) 



where, a(e) and /9(e) = 2arccos(y / e)/7r are functions that 
are equal to unity at e = and are plotted in ref . [2$ ■ The 
time for which the equality holds in J5J will be referred 
to as the speed limit time. 

As shown in ref. Q , it is easy to find models in which 
this limit is approached when one chooses initial states 
in which most of the energy resources are concentrated 
in one subsystem of sufficient simplicity. It also easy 
to construct coherent combinations of basis states that 
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are individually of the form of products of subsystem 
states, in order to achieve fast evolution. But as the 
authors demonstrate, when the system is composed of 
multiple simple subsystems with a Hamiltonian that is 
relatively homogeneous across the range of subsystems, 
and when the system begins in an initial product state, 
then the interaction must generate entanglement among 
the subsystems in order to approach the speed limit for 
the evolution. Ref. Q gives an example in a system of 
N spins, where one of the terms in the Hamiltonian is a 
product of the operators for all spin states, 



h = uj q V(i-^i fe) ) + ^(i-n 



fc=i 



r (*)> 



(3) 



and where the speed limit is realized for choice of pa- 
rameters such that uj >> ^/Nujq. In the present work 
we give some examples of quite different systems that 
can approach the speed limit, again involving quite ho- 
mogeneous Hamiltonians and beginning from a separable 
initial state. But our models are based on a more con- 
ventional pairwise coupling among the spins, 



JV 
i,3 



(4) 



where a+ and <7_ are the usual raising and lowering op- 
erators for an individual spin, and where we will discuss 
different choices for the couplings. 

Within the context of these systems we define an al- 
ternative, and rougher, measure of speed of evolution, 
different from the speed limit rate. This alternative 
definition might, for some purposes, better character- 
ize the system's rate of undergoing macroscopic changes. 
We illustrate by dividing the spin sites into an "up- 
per" tier with j = l,2...,m and a "lower" tier j — 
nl + 1, nl + 2..., nl + n2, where Ui+n2 = N. We consider 
initial states in which all of the upper tier of spins are 
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pointed up and all in the lower tier are pointed down, 

l*(0)> = I T)i ® I T> 2 ® ...| T) ni ® 1 4>n 1+ i ® ...| 4) 

(5) 

Then, to define two different kinds of evolution time, we 
ask the questions: 

1. Following refs 0, 0, for a system with the ini- 
tial state given by JHJ) and governed by one of the 
Hamiltonians that we shall introduce later, what 
is the time required to evolve to a configuration in 
which, 



where 



P(t) = (*(t)\Q 1 \*(t)) <e, 



Q! = |*(0))<*(0)|, 



(6) 



(7) 



and where e is a fairly small number? This is the 
time that defines the speed limit. 

2. For the same system, what is the time required for 
the occupancy of, say, the upper tier of states to 
evolve from the initial condition of all-spins-up to a 
condition in which the upper tier has much different 
occupancy? To quantify this question, we define 



R(t) = (¥(t)|Q 2 |*(t)) 



where 



) 2 = (2n 1 )- 1 ^(l 

i=i 



(8) 



(9) 



and ask, for example, for the time at which R(t) = 
.5. The function R is equal to the average fraction 
of the upper-tier spins that remain pointed up. 

For the Hamiltonian of © in the limit of parameters 
for which the speed limit is approached, the two times 
defined by the above questions are of the same order of 
magnitude. But for our Hamiltonians, they can differ 
greatly, one from the other, in a way that appears to 
depend on the symmetry of the Hamiltonian. 

Before turning to specific cases we note that in the first 
definition of evolution time it would suffice for a single 
spin in the system to be flipped to obtain an orthogonal 
configuration. For our case of a system of a large num- 
ber of spins this would change the average occupancy 
of the upper tier almost not at all. Alternatively, we 
could think of evolution to a product wave-function in 
which each upper-tier single-spin state, originally spin- 
up, is 4 % mixed (in probability) with the spin-down 
state, so that if we had 100 states in the upper tier the 
overlap, P, that enters the first criterion, above, would 
be w exp(— 4); whereas, looking at the second criterion, 



the mixing of the entire upper tier would be only 4%. In 
this case the speed-limit time is much shorter than the 
mixing time. 

We shall be particulary interested in the limit of a large 
number, N, of spins, in our pair-interaction models. The 
speed limit time r will turn out to be proportional to 
iV _1 in these systems, keeping coupling constants fixed. 

MODELS AND THEIR SPEED LIMITS 

In what follows we will consider two possible choices 
for the coupling-constant function in J3J, Xij. 

Case A: Ajj = 1/2. This case has complete permu- 
tational symmetry among the spins. We can write the 
Hamiltonian as, 



H A = ^G(J+J- + J-J+), 



where 



4 = E^ 



(10) 



(ii) 



are the raising and lowering operators for the total spin 
of the system. We specialize to the case ri\ = ri2 = n, for 
simplicity. Then we obtain the values: En = 0, E = ?iG, 
AE = \[2nG . In this case the speed limit will be thus 
given by the term involving E — En in ffi. 



TA > 



7ra(e) 
2nG 



(12) 



Case B: We let the upper set of spins interact with the 
lower set exactly as in case A, but take no interactions 
between pairs of spins both in the upper tier, or both 
in the lower tier. Taking Xij = 1/2 for all upper-lower 
connections and Xij = for the intra-tier connections, 
and defining, 

n 



3=1 
2n 



L ± = £ a 

j=n+l 



(i) 
± > 



(13) 



the Hamiltonian of is, 



H B = G(K + L- + K-L^ 



(14) 



In this case we readily find AE = nG and E=0. Determi- 
nation of the ground state energy En requires a numerical 
calculation which we discuss later. It will turn out that 
in this case the speed limit is determined by the AE term 
in O), 



tb > 



2nG 



(15) 
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Before solving these models we consider some possible 
outcomes for the case of large n. Since G is the only 
dimensional quantity in the theory characteristic times 
are of the form, 



the criterion #1 (i.e. to compare with the speed limit) 
we calculate P(t), defined in ijGj, 



n*)=|E- 



-iGj(j + l)t iy-t-n/2, n/2, 0-|2 
t^n/2, n/2, j J 



(20) 



t = [Gh{n)\ 



(16) 



c; 



mi ,m2 ,m 



for the 



For the evolution to proceed at a finite fraction of the 
speed limits Ijl2|) and Ijl5|) requires h(n) = const x n 
for either case. For small times the simple perturbative 
probability for a single spin to have changed its state is 
of the form, 



prob = const, x G 2 t 2 n 



(17) 



giving a characteristic time for appreciable change of 
single-spin occupation probability of, 



r 2 = const, x (Gy/n)' 



(18) 



In this kind of model, as we indicated earlier, reaching 
a given P(t) < e in the case of large n requires only a 
change of single state probability of order const. /n, and 
the perturbative time to realize this change is, from i|17|) 
of order t\ = const, x (Gn) _1 . Thus we expect to get the 
same dependence on G and n as in the above speed limit 
times, but with a constant that must be calculated. In 
contrast the perturbative estimate of the characteristic 
time for change in R(t) is T2 = const, x (G y/n)^ 1 . When 
we find, for example, that the evolution for large n is 
faster, for example going like 1/n or log(n)/n, we shall 
refer to the phenomenon as a "speed-up" . 



SOLUTION FOR MODEL A 

We take n\ = n% = n to define the initial configuration, 
and we rewrite Ha with J = K + L as, 



H A = G[(K + L) ■ (K + L) — (K 3 + L 3 ) 2 } 



(19) 



Since we are beginning in an eigenstate of (K 3 + i 3 ), 
which is conserved and equal to zero, we can discard the 
final term on the RHS of lfT§|) . Next we note that K ■ K 
and L ■ L are separately conserved under H a ■ The initial 
state is an eigenstate of these two operators with eigen- 
values (rc./2)(n/2+l). The eigenvalues of Ha are given by 
Gj(j + l), associated with the states in which the angular 
momentum K has been added to the angular momentum 
L to give the total angular momentum quantum number, 
j. Thus to follow the time evolution it is only necessary 
to resolve the initial state into eigenstates of J ■ J, and 
let the system develop for some period of time. We then 
take the expectation value of whatever operator, Q, we 
wish to represent the mixing of the system. To address 



where we have used the notation, 
Clebsch-Gordan coefficient. 

For large n and for t < bG^rT 1 ! 2 , where b is some 
number smaller than unity, a sufficiently accurate asymp- 
totic estimate of the Clebsch-Gordan coefficients to be 
used in (|2lH) is given by, 



°n/2, n/2, j 



'(;H-si- 



leading directly to the result, 

P(t) = (1 + G 2 n 2 t 2 ) _1 (l + Q(GV^t) 



(21) 



(22) 



which holds for times less than those of order (Gy^) -1 . 
We have confirmed the limiting result l|22l) with direct 
numerical calculations. In fig.l we plot the function P 
against the dimensionless coordinate, nGt. On the same 
plot we indicate the speed limit rate implied by 112|l . 
plotting the function e(t), where t = wa[e(t)]/2(E — E ), 
E — Eq = nG, and we have used the approximation 
a(e) w /3(e) 2 suggested in ref.Q. We see from fig.l that 
the evolution in the model procedes at roughly half the 
speed-limit rate. 
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FIG. 1: The time evolution of the function P(t) for model A 
(solid line), plotted as a function of the dimensionless variable, 
Gnt. The dashed curve shows the evolution that saturates the 
speed limit. 

We have confirmed the analytic approximations lead- 
ing to (|22|l by direct numerical calculations, and find that 
the curves for P are essentially congruent for all n > 4. 
A numerical calculation over a longer time interval shows 
that after a time of order G" 1 ?^ 1 / 2 , where the approxi- 
mations of (122H fail, the function P(t) remains at a small 
value for an extended time that is of order (Gn) _1 , then 
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reversing its earlier decline in a symmetric fashion; the 
function is periodic with a period 2tt/G, as is clear from 
(HOJl. 

To address criterion #2, based on the evolution of the 
expectation value of the z-component of the total spin of 
the upper tier, we calculate R(t) defined in ©, obtaining, 

*(*) = i-EZ=o£" 

En iGj-(j+l)t^-n/2+fe,ra/2-fe,0f-,-n/2, n/2, 2 
3=0 e n/2, n/2, j n/2, n/2, j \ • 

(23) 

In fig. 2, we show plots of numerical results for n = 
8, 16, 32, 64, 128, 256, in the region in which R(t) declines 
from unity to a value less than .6. 
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FIG. 2: The time evolution of the function R(t) for model A, 
plotted as a function of the dimensionless variable, Gnt, for 
values n = 8, 16, 32, 64, 128, 256. The curves for higher values 
of n lie progressively to the right of the lower n curves. 

The time required to reach R(Gnt) — .6 is strongly 
dependent on n, as we see in fig. 2. If R were to have de- 
veloped at a fixed fraction of the speed limit rate (as did 
P in the previous example), then the value of the dimen- 
sionless variable Gnt required to reach this point would 
have been independent of n . Friedland and Lunardini 
have established that for large n the evolution time 
behaves as (Ci/n) -1 , with a coefficient of order unity; 
so that the characteristic values of Gnt increase as y/n . 
This is essentially the perturbative time dependence that 
we defined above 2 , and is already roughly borne out in 
a fit to the intercepts of the curves shown in fig. 2 with 
the line R = .6. 



2 This result is again obtainable by using the asymptotic forms 
of Clebsch-Gordan coefficients. In place of Friedland and 

Lunardini consider the probability that a single, typical spin is 
flipped, summed over all possible fates for the remaining system, 
and ingeniously reduce the problem in such a way that only a 
single sum analogous to that in 1231 need be performed. 



SOLUTIONS FOR MODEL B 

First we enumerate the set of states that are dynam- 
ically connected, given the initial state JHJ with n% — 
ri2 = n and the Hamiltonian (|14(l . The initial state is 
an eigenstate of K 2 , and L 2 , each of which commutes 
with the Hamiltonian, and each with the same eigen- 
value, (n/2) (n/2 + 1). The total z component "angular 
momentum" , K 3 + L 3 is also conserved and equal to zero 
for our set of states. Any number of the n spins in the 
upper tier, all initially up, may be flipped, leading to n+1 
possibilities for K3 for this tier by itself. For each such 
configuration there is also only a single state of the lower 
tier. Thus we can index the states by the number of flips 
plus one, i, where i takes on the values 1, 2...n + 1. We 
express the operator products that occur in the Hamil- 
tonian in this basis as, 

(i\K-L+\i - 1) = (n - i + l)(i); i = 1, 2...n + 1 
(i - l\K+L-\i) = {n-i + 2)(i - 1); i= l,2...n+ 1 , 

(24) 

which come directly from the standard angular momen- 
tum matrices. We now solve numerically for a (n + 1) 
component wave function *&(t), using the Hamiltonian 
11141 with the substitution (|24|l and the initial condition 
>Fi(0) = 6i t i. We then calculate both of our measures of 
change, 

P(t) = |*r(t)| 2 , (25) 

and 

n+1 

R(t) = Y,\Mt)\ 2 (n-i + l)n- 1 . (26) 
t=i 

We also solve numerically for the ground state energy, 
to confirm that it is the AE condition ofQ that applies. 
We find that the approximation, 

E Q = -OMn - 0.49n 2 , (27) 

fits the energies to within 1% over the range n = 4— 100. 
Thus we have E — Eq > AE over the whole range that 
we consider, and the speed limit is determined from the 
term in involving AE. In fig. 3 we plot P(t) for n = 64 
against the variable gNt. On the same plot we indicate 
the limitation imposed by (|15|l by plotting the function 
e(t) that is the solution to the equation t — n(3(e)/(2 AE) 
where AE = nG. Comparing with the case of model A, 
we see that the evolution proceeds at much more nearly 
the speed limit rate. The plots of P are almost identical 
for all n > 4. 

In fig 4. we show plots of the scaled function R(t) 
for n = 8, 16, 32, 64, 128, 256, 512, again plotting against 
the variable Gnt. Note that the intersections with the 
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FIG. 3: The time evolution of the function P(t) for model B 
(solid line) , plotted as a function of the dimensionless variable, 
Gnt. The dashed curve shows the evolution that saturates the 
speed limit. 

line R = .4 occur at very nearly equally spaced points. 
Equal spacing indicates a mixing time decreasing with 
n as t oc log(n)/?i. This is in contrast to the symmet- 
rical model A in which characteristic time decreases as 
t oc 1/s/n, which we defined as the "perturbative" time. 
Thus we have a large speed-up over the perturbative es- 
timate, as n increases. For the case n = 128, if we de- 
fine the characteristic time according to criterion #2, the 
speed-up over the perturbative rate is a factor of approx- 
imately 4. 
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FIG. 4: The time evolution of the function R(t) for model B, 
plotted as a function of the dimensionless variable, Gnt, for 
values n = 8, 16, 32, 64, 128, 256, 512. The curves for higher 
values of n lie progressively to the right of the lower n curves. 



N , where we taken an initial state which is completely 
separable. We have shown that the time required for 
the squared overlap , P{t), of the evolved state with the 
initial state to become small is indeed proportional to 
TV -1 , with a coefficient which in model A is about twice 
the speed limit value, and in model B is only 10-20% 
greater than the speed limit time, depending on the value 
of e one chooses for purposes of comparison. 

We have also looked at the time scale for change of 
another attribute of these systems, R(t), the fraction of 
upper tier spins that have been flipped. Here the results 
were drastically different for the two models; in model A, 
where there were couplings among all of the spins, the 
evolution time for R(t) is proportional to N~ 1 / 2 , which 
we characterize as the "perturbative" time, rather than 
to N . But in model B, with breaking of permutational 
symmetry 3 the evolution time is reduced to a time of 
order TV -1 log [A], a speed-up to near the A -1 level of 
the speed-limit time. It is tempting to conjecture that 
with more breaking of the symmetry we could lose the 
log, and perhaps even approach speed-limit rates. The 
curves shown in ref. ||, however, for the case of com- 
pletely broken symmetry, were inconclusive in this re- 
gard, due to computational limitations. More generally, 
it would seem to us to be of some interest first to define 
"macroscopic change" more precisely in a large N limit, 
and then to try to find general results, analogous to the 
speed-limit results, on the limitations for the rates for 
this macroscopic change. 

We should note that we have tested models in which 
n\ 7^ n- 2 and find results qualitatively the same as those 
reported for the two models treated in this paper. We 
set out, in calculating R(t) to get an estimate of the 
rate at which the average total z component spin of the 
upper tier, A3, changes. But see that in model B the 
measure R(t) is very nearly zero at the first minimum, 
indicating that almost all the spins have been flipped. 
We could compare the results, where (up to logarithms) 
the characteristic time for evolution is of order (GA) -1 , 
to the results for the simplest Hamiltonian that will flip 
each spin in time (GA) _1 , namely, 

Hc = \GNj^o^. (28) 

i=l 

Measuring the "energy resources" of this Hamiltonian by 
AE = (tt/2)GA 3 / 2 , and comparing to AE = G(A/2) for 



DISCUSSION 

We have considered large spin systems with couplings 
among pairs such that the quantum speed limit time is 
proportional to the inverse of the number of spins, r oc 



3 Note that we broke the spin states into two sets both for the 
purpose of breaking the permutation symmetry by setting to 
zero the intraset interactions, and for the purpose defining which 
states were initially up and which down. We used the same 
division for the two purposes only to be able (for large N) to 
reduce the evolution problem to calculable dimensions. This is an 
example of how very minimal is our breaking of the symmetries 
of the symmetric model. 
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the Hamiltonian {OJ, we conclude that the model with 
pair couplings and consequent many-body entanglements 
does indeed succeed in rotating the complete ensemble of 
spins in a very efficient way, for large n. 

Finally, we comment on the question as to how much 
substance is contained in the statement that "entangle- 
ment is the key to speed-up of rates" a paraphrase of a 
point emphasized in and also in a quite different con- 
text in [Sj. Of course, states do become entangled during 
the evolution under the interactions that we have taken 
in our models. But there is no general agreement known 
to us on the quantification of entanglement in the case of 
multicomponent systems, even in the case of pure states. 
We can, however, perhaps rather arbitrarily, consider the 
entanglement of two systems, U and L, where U contains 
all of the upper tier spins and L contains the lower tier 
states. Then a measure of entanglement between the two 
systems is given by the entropy of entanglement obtained 
by tracing out one subsystem's coordinates in the density 
matrix, say the U subsystem, to define a reduced density 
matrix pl for the L subsystem, and then calculating the 
von-Neumann entropy corresponding to Pl |l0| ] [TlT ] , 

S e = -Trp L log[p L ] . (29) 

We have calculated S e as a function of time for models 
A and B, using the wave- function solutions underlying 
the calculations discussed above. At a time Gnt = 1, at 
which P .5 in either model, we find that, so defined, the 
entanglement is about 30% more for the case with broken 
symmetry than for the completely symmetric case, and 
relatively independent of the number of particles. Thus 
in our models there is an perhaps an indication of a cor- 
relation between the rate of growth of entanglement and 
the efficiency of transformation. Of course, in the pure 
states that we have considered, the entanglement, as de- 
fined above, can both increase and decrease in time, and 
comparisons will depend on the time of the sampling. 
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